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Abstract. Although soft X-ray absorption features in 
AGN were discovered almost ten years ago, the nature and 
location of the gas creating them has remained controver- 
sial. However, by making use of the newly-available high- 
resolution spectra provided by XMM-Newton and Chan- 
dra, we should be able to make substantial advances. The 
first such spectra indicate that multiple ionization states 
coexist in the absorber; this is a natural consequence of 
photoionization physics. Photoionized evaporation in the 
presence of a copious mass source locks the ratio of ioniz- 
ing intensity to pressure to a critical value. A broad range 
of temperatures can all coexist in equilibrium for this value 
of the ratio of ionizing intensity to pressure. Consequently, 
the flow is expected to be strongly inhomogeneous in tem- 
perature. The inferred distance of this material from the 
source of ionizing radiation depends on how much mat- 
ter exists at the highest-obtainable temperature. This dis- 
tance can be measured by monitoring how ionic column 
densities respond to changes in the ionizing continuum on 
timescales of days to years. 



1. The Unanswered Questions about Warm 
Absorbers: How Chandra and XMM Can Help 

Although absorption features due to highly-ionized species 
were first clearly seen in AGN spectra taken with ROSAT 
almost ten years ago (Turner et al. 1993), most of the 
fundamental facts about them remain obscure today. This 
is true at even the most basic, empirical level. 

1.1. The ionization distribution 

For example, we do not even know which ionization states 
are present in the absorbing gas. Because increasing levels 
of ionization tend to move atomic features higher in energy 
and then remove them altogether, the character of ioniza- 
tion states possible to see in X-ray spectra is strongly in- 
fluenced by the specifics of instrumentation. Consequently, 
our efforts to determine how much matter there is in which 
ionization state suffer from extremely strong instrumental 
selection effects. 



In ROSAT-data, the spectral resolution and sig- 
nal/noise were so poor that the only feature strong enough 
to be seen was the O K-edge, and the ionization state of 

could be determined only roughly. Labelling ionization 
level by the ionization parameter £ = Li on /nr 2 (where 
Lion is conventionally defined as the luminosity between 

1 and 1000 Ryd, n is the H-nucleus density, and r is the 
distance to the source of radiation), this sort of data was 
capable only of identifying the presence of matter with 
£ < 30, and perhaps distinguishing ionization parameters 
at the highest end of this range from those an order of 
magnitude or more smaller. 

Unfortunately, ROSAT was was totally incapable of 
detecting any matter more strongly ionized because higer 
ionization destroys OVIII. Other elements' K-edges could 
not be seen in ROSAT data either because their energy 
was outside its limited range (e.g., Fe) or because their 
abundance was too small to permit detection with such 
crude data. To put this (semi)-quantitatively, in gas with 
solar abundances most of whose O has not been stripped, 
O K-edges were detectable by ROSAT when the column 
density N ^ 10 21 cm" 2 . Even if the ROSAT energy 
range extended to the K-edge of ionized Fe at 7-9 keV, 
a column density almost three orders of magnitude larger 
would have been necessary in order for a significant edge 
to have been formed. In this sense, ROSAT-level data was 
roughly 1000 times less sensitive to highly-ionized than to 
weakly-ionized gas. 

Data of the sort produced by ASCA was better, but 
not greatly. Although its resolution and signal/noise were 
good enough to permit estimates of the depths of the OVII 
and OVIII edges separately (e.g., Reynolds et al. 1995), it 
could see no other edge features reliably. The next-deepest 
K-edges (Si, S, Ar, etc.) are so much smaller in opacity 
that they would be detectable in ASCA data only for much 
greater absorbing columns. In principle, ASCA's nominal 
energy range included the Fe K-edge, but the effective area 
there was so small as to make searching for this edge im- 
practical. Thus, like ROSAT, ASCA could detect matter 
with £ <^ 30, but could improve upon ROSAT by permit- 
ting £ ~ 30 to be distinguished from, say, £ = 3. Again like 
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ROSAT, it was effectively blind to more highly-ionized 
matter. 

The advent of the grating spectrometers on Chandra 
and XMM has already dramatically changed this frustrat- 
ing situation, and will undoubtedly improve matters fur- 
ther in the near future. Their tremendously finer spectral 
resolution permits study of line transitions, not just edges, 
while their much greater signal/noise allows one to search 
for much weaker features. With data of this quality, it is 
possible to search for species expected to be common up 
to £ ~ 1000 even when the column density of that gas is as 
small as <; 10 21 cm~ 2 — and we are starting to find them. 
Several different recent warm absorber studies employing 
grating data report that more than one ionization compo- 
nent is detected, with £ ranging anywhere from a few tens 
to - 1000 (Kaastra et al. 2000, Collinge et al. 2001, Kaspi 
et al. 2001). 

Thus, in the next few years we may reasonably hope 
to make much progress on the most basic empirical ques- 
tion of warm absorber studies: measuring the distribution 
dN/dt 

1.2. Location and origin 

After identifying the absorber's ionization state, the next 
obvious question is where to place it in the AGN system. 
Proposed distances have ranged all the way from inside 
or near the broad-line region at 0.01 - 0.1 pc from the 
nucleus (Reynolds et al. 1995, George et al. 1998, Netzer 
et al. 2001) to many tens of pc away in the narrow-line re- 
gion (Bottorff et al. 2000). Some have argued that there is 
absorbing material across this entire range of radii (Otani 
et al. 1996, Morales et al. 2000). Others (e.g., Krolik & 
Kriss 1995) have argued that the absorbing gas should 
be identified with the scattering region posited by Seyfert 
galaxy unification schemes, perhaps ~ 1 pc from the nu- 
cleus when L ion ~ 10 44 erg s _1 . 

The state of affairs regarding the absorbing matter's 
origin is even worse. A few of the diverse suggestions that 
have been floated include: evaporation off "bloated stars" 
in the broad-line region (Netzer 1996); gas evaporated off 
the torus obscuring the nucleus (Krolik & Kriss 1995); and 
a wind driven off the accretion disk (Elvis 2000; Bottorff 
et al. 2000). The absorbing gas's ultimate fate has been 
left wholly unconsidered. 

This, too, is a topic on which we may hope for en- 
lightenment from new spectroscopy data. As has been re- 
marked in numerous papers (Reynolds et al. 1995, Krolik 
& Kriss 2001, Netzer et al. 2001), variability in absorp- 
tion features can strongly constrain the location of the 
gas responsible. The virtue of the new data sets is that, 
by making so many features observable, our odds of seeing 
variability in one or more of them are greatly improved. 
Moreover, because the new data are less subject to ion- 
ization selection effects (as discussed in the previous sub- 



section), the location bounds placed by variability can be 
applied more generally. 

1.3. Kinematics 

Independent of where the absorbing gas is located, we'd 
also like to know which way it's moving. For this issue, 
high-resolution spectroscopy is essential. If the absorbing 
gas possesses ions with ionization potentials ~ 1 keV (as is 
necessary to have numerous soft X-ray features), its tem- 
perature cannot be much greater than ~ 10 7 K. Transonic 
motions at this temperature would then be ~ 300 km s _1 
or slower. It immediately follows that spectral resolution 
AE/E 1000 is a prerequisite for beginning to detect 
motions at this speed. 

Because resolution of this caliber is exactly what Chan- 
dra and XMM have now delivered, we are beginning to 
learn about this question. Kaspi et al. (2001), for exam- 
ple, have shown that the absorber in NGC 3783 is moving 
toward us (relative to the nucleus) at ~ 500 km s _1 . 

1.4- Geometry 

Although the depth of an edge may tell us the total col- 
umn density of ions along the line of sight, it does not tell 
us whether they are all clumped together, or spread out 
smoothly, or clustered at a variety of locations. And, of 
course, absorption tells us nothing about the existence of 
matter along other directions. 

Other constraints are necessary to make progress an- 
swering questions of this sort. For example, resonance line 
emission components offset toward the red from absorp- 
tion components are the classic signature of outflows. The 
relative equivalent width of the emission and absorption 
components reveals the ratio between the mean optical 
depth on lines of sight other than ours to the ratio on 
the line toward us. Similarly, the relative velocity widths 
further constrain the gas's geometric disposition. Making 
use of the emission/absorption equivalent width indicator, 
Kaspi et al. (2001) have already suggested that the absorb- 
ing gas in NGC 3783, while occupying a large part of solid 
angle around the nucleus, does not completely cover it. 

dumpiness in the gas can be analyzed in any of sev- 
eral ways, given quality spectroscopy: Emission lines that 
are either collisionally-excited or created by recombina- 
tion yield measurements of the gas's emission measure; 
combined with independent density estimates, one can 
then infer the total emitting volume and compare it to 
the available volume in the region. Alternatively, spectro- 
scopic data that requires gas at several different ionization 
levels co-existing in more or less the same location would 
strongly suggest a clumpy, inhomogeneous environment. 
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1.5. Dynamics 

Given our ignorance about the warm absorber's kinemat- 
ics, location, and origin, it should come as no surprise that 
almost nothing is known with confidence about its dy- 
namics. We have as yet advanced hardly at all past listing 
the "usual suspects": gravity, thermal presure gradients 
(Balsara & Krolik 1993), radiation pressure (Chelouche & 
Netzer 2001, Morales & Fabian 2002), and magnetic fields 
(Bottorffet al. 2000). 

2. Warm Photoionized Equilibria 

To make further progress, it is helpful to consider a few 
constraints imposed by the basic facts of photoionization 
physics. We begin with the most elementary of these, a 
relationship that links the ionization parameter, the dis- 
tance, the gas column density, and its volume filling factor. 
From the definition of £ and the fact that the column den- 
sity N — nAr, where Ar is the total distance along the 
line of sight occupied by the absorbing gas, we immedi- 
ately see (as shown by Turner et al. 1993) that 



Ar 



r = 



a) 



Because Ar/r cannot possibly be greater than unity, 
the inferred ionization parameter and column density of 
an absorber, when combined with the luminosity of the 
ionizing source, place an upper bound on its distance from 
that source: 

T ^ Tmax — AJ -j. • (2) 



For typical numbers, r max ~ 30L io „ i 447V^ 1 (^/100) _1 pc. 
Here we have scaled the ionizing luminosity to 10 44 erg s _1 
and the column density to 10 cm . For fixed r ma;r , the 
more highly clumped the gas is, the smaller the radius at 
which it might be found. 

A second constraint is provided by the nature of 
photoionization equilibrium. If the gas's cooling time is 
shorter than the time required for it flow through the ab- 
sorbing region, we can expect that it is in energy balance 
with the radiation. If that is the case, its temperature is 
a unique function of its ionization parameter £. 

If we wish to link photoionization physics with fluid 
dynamics, it is more convenient to write the ionization 
parameter in a related, but significantly different form: 
S = ^/(i-KckT). The factor of At: ck makes S renders it 
dimensionless; more significantly, by dividing by the tem- 
perature, a given S at any particular radius corresponds to 
a certain pressure. Viewed from the opposite perspective, 
if the boundary condition is a fixed pressure and the gas is 
in energy balance, 5 is determined, but T is uniquely de- 
termined only when T(S) is single- valued. Although T(£) 
is, in general, unique, when the ionizing continuum has 
a form typical of AGN, T(S), in general, is multi-valued 
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SXVI, ArXVIII, FeXXIV, FeXXV 

NeX, MgXII, SiXIV, ArXVII, 
FeXIX, FeXX 



OVIII, NeX, MgXI, MgXII, 
SiXIII, FeXVII 



OVII, OVIII, NelX 



log(E) 

Fig. 1. The temperature equilibrium curve for two guesses 
(Krolik & Kriss 2001-solid curve, Kaspi et al. 2001-dashcd 
curve) about the unknown shape of the EUV spectrum in 
Seyfert galaxies. Although the value of S at which the 
nearly vertical rise between 5 x 10 4 K and 1 x 10 6 K takes 
place depends on the shape of the spectrum, the character 
of the rise is a persistent feature of all reasonable AGN 
spectra. The location of several values of £ is shown, as 
are a few of the dominant ionization stages associated with 
each of the marked values. 



in the range of temperatures between ~ 5 x 10 4 K and 
£ 10 6 K (Krolik, McKee & Tarter 1981; fig. 1). As a re- 
sult, the boundary conditions may not suffice to deter- 
mine the temperature; its history (and possibly the influ- 
ence of thermal instabilities) is instead the final deciding 
factor. The fact that the thermal balance may be multi- 
valued is central to the study of warm absorbers because 
the multiple- valued solutions occur in exactly the range of 
conditions capable of making soft X-ray features. 

When a detailed calculation is made of the thermal 
equilibrium curve, it is necessary to make several guesses 
about matters that we don't truly know and yet have no- 
ticeable effect on the result. The distribution of elemental 
abundances is the first of these. Generally-speaking, solar 
abundances are assumed; other choices can alter the de- 
tailed shape of the curve, but within quite wide bounds 
have at most modest effect on the shape of T(E). More 
important to the thermal balance curve is the shape of 
the ionizing continuum. In low-redshift objects, we can 
observe the continuum at energies just below 1 Rydberg 
and again above ~ 500 eV, but the most important seg- 
ment of the ionizing continuum is precisely in the band 
from 13 cV to ~ 500 eV that we can never directly see. 
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The best we can do is smoothly interpolate from either 
side. Fortunately, as shown in Figure 1, different plausi- 
ble guesses do not result in qualitatively different results. 
Both equilibrium calculations shown in the figure were 
based on continuum shapes meant to apply to the Seyfert 
galaxy NGC 3783 (one in Kaspi et al. 200f , the other in 
Krolik & Kriss 200 f ). Although rather different, both pre- 
dict the same general sort of behavior — a nearly vertical 
rise in temperature from ~ 5 x f0 4 K to almost fO 6 K. 
The value of E where this takes place is called E c . 

In fact, it is worthwhile to focus greater attention on 
the "vertical branch" of the photoionized thermal equi- 
librium curve because it is significant for several reasons. 
The first, as just mentioned, is that the range of £ inferred 
from the ionization stages seen in warm absorber spectra 
is ~ 10 - 1000; this is almost exactly the range traversed 
on the vertical branch. In other words, if warm absorbers 
are in a state of photoionization equilibrium and thermal 
balance, they "live" on the vertical branch. 

Second, the vertical branch has a rather special ther- 
modynamic property: it is marginally stable to isobaric 
thermal perturbations. Were dT/dE > 0, constant pres- 
sure perturbations would be stable because increasing the 
temperature from an initial point on the equilibrium curve 
must result in net cooling; conversely, were dT/dE < 0, 
the same sort of perturbation grows exponentially be- 
cause the excursion causes net heating. However, on the 
vertical branch dE/dT = 0: in this case, if the tempera- 
ture is perturbed while the pressure is fixed, the gas re- 
mains in a state of zero net radiative energy gain or loss. 
Consequently, there is no temperature "feedback", either 
positive or negative. Perturbations are neither strongly 
damped nor strongly amplified — they simply do whatever 
the outside agency creating them demands. Although it is 
true that the curve of thermal balance can never be pre- 
cisely vertical, to the extent that it is nearly so, the net 
cooling or heating that is engendered by the perturbations 
is small, and the magnitude of growth permitted in a re- 
gion of weak instability is capped at a low level as a new 
segment of the equilibrium curve is quickly reached. 

Third, far from being a special case, the value of 5 
at which the vertical branch occurs is actually a pre- 
ferred value of the pressure-based ionization parameter. 
Whenever large amounts of cool material are exposed to 
a strong ionizing flux (as might readibly be imagined near 
an AGN), the low-density portion of this cool matter finds 
itself rapidly ionized and heated. If the heating time is 
short compared to the flow time, the ambient pressure 
rises equally rapidly. At fixed distance from the source of 
ionizing radiation, this means the value of E surrounding 
the cool matter rapidly falls. Nothing stops this fall in E 
until it reaches E c , where the net energy exchange via ra- 
diation switches from heating to cooling. Similarly, if the 
pressure overshoots to the net-cooling side of the vertical 
branch, matter cools rapidly, losing pressure, until balance 
is achieved — at exactly E c . 



This portion of the argument may be summed up very 
simply: Given that, by definition, warm absorbers are ex- 
posed to AGN fluxes, it makes sense to suppose that they 
are photoionized. Let us suppose for the moment that the 
gas is in photoionization equilibrium and radiative bal- 
ance. Then the particular ions observed (OVII,VIII; NeX; 
MgXII; SiXIV; etc.) indicate that the gas lies on the verti- 
cal branch of the T(E) curve. Because that branch allows a 
wide range of temperatures to coexist at a single pressure, 
it follows that the material is likely to contain multiple 
sub-regions whose temperature could be anywhere within 
the range found on the vertical branch of the equilibrium 
curve. Moreover, pressure balance tends to force photoion- 
ized gases in the presence of a copious mass source toward 
E c . That is, it is a natural consequence of photoionization 
physics that we should see gas with E = E c and that it 
should be inhomogeneous in temperature. 

3. Global Picture 

These special thermodynamic properties of the vertical 
branch find a natural application in the context of AGN 
warm absorbers. These are frequently seen in type 1 
Seyfert galaxies, a variety of AGN in which we have ex- 
cellent reason to believe that extremely optically thick 
obscuring matter wraps around the nucleus with roughly 
toroidal geometry (Antonucci 1993, Krolik 1999). Where 
the ionizing radiation of the nucleus strikes the inner edge 
of the obscuring stuff, exactly the process described in 
the previous paragraph can be expected to take place: co- 
pious amounts of matter should be ionized and heated, 
maintaining 5 ~ E c . Until the temperature rises quite 
high, the heating time is far shorter than the flow time: as 
shown in detail in Krolik & Kriss (2001), if the flow speed 
is of the magnitude observed (~ 500 km s _1 ), this crite- 
rion is met for T < 4 x lO 6 ^ 1 ^ 2 K, for r pc the distance 
from the central source in parsecs. 

As the newly-liberated material is heated, it will ex- 
pand in order to maintain pressure balance. One might 
then ask, "Will the expanding parcels eventually be con- 
fined? Or will they spread out and fill the entire volume 
of the region?" Suppose that the matter we see by its 
"warm absorber" features is clumpy. If the volume out- 
side the clumps has lower pressure, the clumps will ex- 
pand at roughly their internal sound speeds. The ratio in 
column density between the low-pressure background and 
the high-pressure absorbing lumps may then be estimated 
as 

N bkgd /N abs ~ (r/Ar abs )(c s /v flow ), (3) 

where the absorber has radial thickness Ar abs and sound 
speed c s . Because the sound speed at ~ 10 5 K (the typ- 
ical temperature inferred for warm absorber regions) is 
~ 50 km s _1 , it is clear that the column density of the 
background is automatically close to the column density 
of the absorber if there is strong clumping. That is, one 
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must expect the background to be a significant part of the 
absorbing structure. Put yet another way, it is probably 
best to consider the entire system as volume-filling, albeit 
one that is likely to be strongly inhomogencous. 

Although the ionization parameter S is tightly con- 
strained by photoionized thermodynamics, the tempera- 
ture, at least within the range permitted at S ~ S c , is 
hardly constrained at all. We should expect, then, that 
the evaporating matter in this region might exist any- 
where in the temperature range 5 x 10 4 - ~ 10 6 K, i.e., 
30 <; £ <; 1000. The issue is not "What is the favored 
value of T or £?" , but "What is the distribution of T and 
£ within this broad range?" Small accidents having to do 
with irregularities in the source of matter (the obscuring 
torus?) or "bumps" in the outflow of heated gas could 
affect this distribution, so different objects could easily 
differ in detail. 

At fixed pressure, the hottest regions have the least 
density and so, per unit mass, occupy the largest volume. 
For this reason, one might expect most of the volume to be 
occupied by relatively hot gas, but it is likely studded by 
smaller regions of higher density and lower temperature. 
Each of these sub-regions can easily change its tempera- 
ture over time. 

4. Implications 

This conceptual picture has numerous implications for the 
character of warm absorbers. The most basic, of course, is 
that the features we see should originate in ions found in 
photoionized equilibrium over a span of a factor of 20 or 
so in temperature. There are also other, more quantitative 
implications. 

4-.1. Location 

Returning to equations 1 and 2, we can rewrite them as 
r (£ ) 1 

1 max \smax J 



where ^ max is the highest £ (corresponding to the highest 
T) attainable in a flow time. Because the gas at £, m ax is 
the hottest that we can expect to find in the absorber 
region, it should occupy the largest part of the volume. 
That is, for the hottest gas Ar/r ~ 1, while cooler gas is 
strongly clumped. Scaled in terms of typical numbers, the 
predicted distance 



n Li on 44 



(5) 



where £, ma ,x has been set at 1000. 

A few parsecs is, of course, the scale of the inner torus 
(as seen, for example, in near-infrared imaging — Thatte 
et al. 1997, Marco & Alloin 2000 — or H 2 maser spots— 
Greenhill et al. 1996). It is likewise the expected scale of 
the polarizing reflection region that permits us to see type 
1 Seyfert nuclei in type 2 Seyfert galaxies. 



4.2. Dynamics 

With an estimate of the distance to the nucleus, the 
strength of gravity may be more quantitatively estimated. 
Making use of equation 5, we find that the free-fall speed in 
the absorber region is ~ 30(L/L E y 1/2 N^ 2 km s" 1 . Here 
Le is the Eddington luminosity and we have supposed 
that most of the bolomctric luminosity is in the ionizing 
band. Because the sound speed in the hottest phase of the 
absorber is ~ 100 km s _1 , we can expect that thermal 
pressures alone may be able to drive an outflow unless 
L/Le <C 1. Whether radiation or magnetic forces sup- 
plement the thermal pressure gradient remains an open 
question, but the fact that we see blue-shifted absorption 
lines indicating an outflow should now be no surprise. 

4-3. Geometry 

If the warm absorber is formed by this sort of photoionized 
evaporation off the inner edge of the obscuring torus, the 
torus itself blocks the warm gas from expanding to the 
side. A biconical expansion and outflow results. For this 
reason, we might expect the warm absorber to occupy fully 
only a part of solid angle around the nucleus. Moreover, 
depending on our viewing angle, parts of the obscuring 
matter might block our view of selected portions of the far 
side of the outflow. Thus, emission line equivalent widths 
should indicate substantial, but incomplete, solid angle 
coverage — as is tentatively seen in NGC 3783 (Kaspi et 
al. 2001). 

4-4- Variability and location 

This picture also makes specific predictions about vari- 
ability in absorption features. The column densities we 
see are, of course, controlled by ionization balance in the 
line-of-sight gas. The governing equation can be written 
(in abridged form) as 



din rii 
dt 



m 



n e cx rpr — — h • • ■ 



Fx® ion W>e®-rec 



(6) 
(7) 



where the density of ionization stage i is m, i ± 1 corre- 
spond to the stages once more or less ionized, F x is the 
photon flux at the relevant ionization edge, 0~i O n is the 
appropriately averaged ionization cross section, and the 
primes and double primes denote the quantities relevant 
to the i + 1 and i — 1 ionization stages. When n i+1 ;§> 
(as is often true for the particular ionization stages seen), 
two terms dominate the rest, so that 

d\nni .. 

d - n e a rec n i+1 /n t - F x a ion . (8) 

The characteristic time for effecting changes in the ab- 
sorbing column is then 

„2 



ti< 



10 b 



-'ion, 44 



(9) 
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Different species have varying ionization timescales, de- 
pending on the details of their atomic structures. For 
example, if one examines the more abundant species in 
the relevant ionization conditions, SiXIV, FeXXV, and 
FeXXVI tend to have longer ionization times, while OVI 
and the ions around FeXX tend to have somewhat shorter 
timescales. The range is not terribly large, however: the 
±1 in the exponent in equation 9 more or less comprises 
the range for the species of greatest interest (Krolik & 
Kriss 2001). 

Timescales of days to years are, of course, very conve- 
nient timescales for human observations. Because we know 
which species should vary on which timescales (parame- 
terized by the ionizing flux in the absorbing region), mon- 
itoring spectroscopy should be a very powerful tool for 
constraining the location and conditions of the absorbing 
gas. 
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4-5. How would the warm absorber look when seen from 
the side? 

Finally, it is striking how much the warm absorber region 
resembles the reflection region that allows us to see the 
nuclei of type 2 Seyfert galaxies. Both regions appear to 
lie at distances ~ pc; both are highly-ionized; both are 
flowing outward at speeds of several hundred km s _1 ; the 
total column densities in both are ~ 10 22 - 10 23 cm~ 2 . It 
would appear to make a great deal of sense to suppose that 
they are the same structure, simply viewed at different 
angles. 
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